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ABSTRACT 

We present millimeter and radio observations of 13 SDSS quasars at reshifts 
z~6. We observed eleven of them with the Max-Planck Millimeter Bolometer 
Array (MAMBO-2) at the IRAM 30m-telescope at 250 GHz and all of them with 
the Very Large Array (VLA) at 1.4 GHz. Four sources are detected by MAMBO- 
2 and six are detected by the VLA at > 3a level. These sources, together with 
another 6 published in previous papers, yield a submillimeter/millimeter and 
radio observed SDSS quasar sample at z~6. We use this sample to investigate 
the far-infrared (FIR) and radio properties of optically bright quasars in the early 
universe. We compare this sample to lower redshift samples of quasars observed 
in the submillimeter and millimeter wavelengths ((sub) mm), and find that the 
distribution of the FIR to B band optical luminosity ratio (Lfir/Lb) is similar 
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from z~2 to 6. We find a weak correlation between the FIR luminosity (Lpm) 
and B band optical luminosity (Lb) by including the (sub)mm observed samples 
at all redshifts. Some strong (sub) mm detections in the z~6 sample have radio- 
to-FIR ratios within the range defined by star forming galaxies, which suggests 
possible co-eval star forming activity with the powerful AGN in these sources. We 
calculate the rest frame radio to optical ratios ( 4 = L,^^i.4ghz/L ^^^q^) for all 
of the VLA observed sources in the z~6 quasar sample. Only one radio detection 
in this sample, J083643. 85+005453. 3, has R* 4 ~ 40 and can be considered radio 
loud. There are no strong radio sources (R^ 4 > 100) among these SDSS quasars 
at z~6. These data are consistent with, although do not set strong constraints 
on, a decreasing radio-loud quasar fraction with increasing redshift. 

Subject headings: galaxies: quasars — infrared: galaxies — radio continuum: 
galaxies — galaxies: starburst — galaxies: high-redshift 



1. Introduction 

The tight correlation between the mass of supermassive black holes (SMBH) in the 
centers of galaxies and the bulge mass/velocity dispersion is well documented in the local 
universe (Magorrian et al. 1998; Marconi & Hunt 2003; Ferrarese & Mcrritt 2000; Gcbhardt 
et al. 2000; Tremaine et al. 2002) and suggests that the evolution of SMBH and spheroidal 
galaxies is coupled even at high redshift. Due to the negative K-correction, submillimeter 
and millimeter ((sub) mm) observations become an efficient way to study the host galaxies of 
high redshift quasars (Blain & Longair 1993). Such observations probe the rest frame FIR 
emission from the dust components in the interstellar medium (ISM) of these objects, and 
thus provide a unique chance to test co-eval black hole and bulge formation by searching for 
massive starbursts in the host galaxies of high redshift optically bright quasars. 

The (sub) mm and radio properties of some optically bright quasars at high redshift 
are discussed in a number of papers. Observations of large samples of optically selected 
quasars from z~1.5 to 5 (Omont 2001, 2003; CariUi et al. 2001; Isaak et al. 2002; Priddey 
et al. 2003a; Beelen et al. 2004) result in a (sub) mm detection rate of about 30% at 
mJy sensitivity. The derived far-infrared (FIR) luminosities of the (sub) mm detections arc 
typically ~ 10^^ Lq and imply dust masses > 10® Mq. Comparing samples observed with the 
Max-Planck MiUimeter Bolometer (MAMBO) at the 30-meter IRAM telescope at redshift 2 
and 4, Omont et al. (2003) found that the FIR luminosities of the optically bright quasars 
in these samples do not evolve with redshift. Statistical tests for these (sub) mm observed 
quasars show a weak correlation between the optical and FIR luminosities (Omont et al. 
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2003; Beelen 2004; Cox et al. 2005), which is argued as an evidence for FIR emission from 
warm dust heated by star formation. The inferred star formation rate is ~ 10^ Mq yr~^ 
when assuming a normal initial mass function (IMF). This interpretation is supported by 
the fact that the radio-to-FIR ratios for the FIR-luminous sources are consistent with the 
range spanned by star forming galaxies (Carilli et al. 2001; Petric et al. 2006). 

The Sloan Digital Sky Survey (SDSS, York et al. 2000) has identified 19 bright quasars 
at z~6 (Fan et al. 2000, 2001, 2003, 2004, 2006a). These quasars are imique, in that they 
are undergoing rapid accretion onto SMBHs with masses > 10^ Mo, within 1 Gyr of the 
Big Bang, an epoch approaching cosmic reionization (Fan et al. 2006c). These Z'^6 quasars 
have been observed at all wavelengths from the X-ray to the radio (eg. Bechtold et al. 
2003; Pentericci et al. 2003; Shemmer et al. 2006; Jiang et al. 2006a; Petric et al. 2003). 
The X-ray to near-infrared (NIR) observations imply Spectral Energy distributions (SEDs) 
similar to those of local quasars. The integrated bolometric luminosities exceed > 10^^ L©, 
and the black hole masses are > 10^ M© (Willott et al. 2003; Iwamuro et al. 2004; Jiang 
et al. 2006a). Seven z~6 SDSS quasars (Fan et al. 2000, 2001, 2003) were observed at 350 
GHz with the SCUBA camera at the James Clerk Maxwell Telescope or at 250 GHz with 
MAMBO. Three are detected by SCUBA (Priddey et al. 2003b; Robson et al. 2004) and 
two are detected by MAMBO (Bertoldi et al. 2003a). Deep VLA observations at 1.4 GHz 
were made for six sources, resulting in detections of two of them (Petric et al. 2003; Carilli 
et al. 2004). 

SDSS J114816. 64+525150. 3 at ^ = 6.42 is the best studied FIR luminous source at 
z~6. It was detected with SHARC II, the bolometer camera at the Caltech Submillimeter 
Observatory at 350 /im, SCUBA at 450 /im and 850 /xm, MAMBO at 1.2 mm, and the VLA 
at 1.4GHz (Beelen et al. 2006; Robson et al. 2004; Bertoldi et al. 2003a; CariUi et al. 2004). 
Fits to the SED from the FIR to radio imply a warm dust component with a temperature 
of 55K (Beelen et al. 2006). The corresponding FIR luminosity is 2.2 x 10^^ with an 
inferred dust mass of 4.5 x 10^ M©. It is also the only z~6 source that has been detected in 
the CO and [C II] 158 fim emission lines (Bertoldi et al. 2003b; Waher et al. 2003; Walter 
et al. 2004; Maiolino et al. 2005). These CO observations reveal a large mass of molecular 
gas (~ 2 X 10-*^° Mq) in the host galaxy. 

The (sub)mm through radio observations of J114816. 64+525150. 3 argue for a massive 
starburst in its host galaxy, because (i) the FIR to radio ratio of this source follows that of 
typical star forming galaxies in the local universe (Yun et al. 2001; Carilli et al. 2004), (ii) 
the huge amount of molecular gas in its host galaxy can provide the required fuel for massive 
star formation, and (iii) the star formation rates indicated by both FIR luminosity and [C 
II] 158//m fine luminosity are ~ 10^ Mq yr'^ (Bertoldi et al. 2003a; Maiolino et al. 2005). 
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Nineteen z~6 SDSS quasars have been published to date (Fan et al. 2006b). These are 
an optically selected sample at the highest redshift. We are pursuing a series of (sub)mm 
through radio studies on this sample in order to (i) find the general dust and gas properties in 
the host galaxies of these highest redshift quasars, and (ii) search for star formation activity 
co-eval with the rapid growth of SMBHs in the early universe. 

In this paper, we present new MAMBO-2 250 GHz observations of eleven, and VLA 
1.4GHz observations of thirteen, z~6 SDSS quasars. Then, together with previously pub- 
lished results, we discuss the general far-infrared and radio properties of the optically selected 
quasars at z~6, focusing on luminosity correlations and evolution. We will also discuss the 
radio loud fraction of quasars. We will present further analysis of the FIR to radio SEDs 
and discuss possible star forming activity in a second paper (paper II). The sample and 
observations are described in section 2, and results are presented in section 3. In section 4 
and 5, we analyse and discuss the general properties of FIR and radio luminosities for the 
entire sample and give the conclusion in section 6. We adopt a concordance cosmology with 
Hq = 71 kms~^ Mpc~^, Qm = 0.27 and Qx = 0.73 throughout this paper. 



2. Sample and Observation 

The sample of z~6 quasars are selected from about 6600 deg^ of imaging data from the 
Sloan Digital Sky Survey (Fan et al. 2006b). Fan et al. (2000, 2001, 2003, 2004, 2006a) 
selected sources with very red i — z colors as z> 5.7 quasar candidates. These candidate 
sources were followed up with high quality spectroscopy. Nineteen z~6 quasars have been 
published by the SDSS survey (Fan et al. 2000, 2001, 2003, 2004, 2006a) with redshifts 
ranging from 5.74 to 6.42 and rest frame B band optical luminositie^ from 10^^'^ Lq to 
IO^'^ '^Lq. Most of the 13 objects we observed in this work were discovered in the past two 
years (Fan et al. 2004; 2006a) and have not been observed in the (sub)mm or radio bands 
before. 

We compare the (sub)mm properties of the z~6 objects with luminous {Lb > 10^^'^ Lq) 
quasars at lower redshifts from the literature. We define a low redshift group, 1.5 < z < 3.0, 
using a sample from Omont et al. (2003) containing radio-quiet quasars with B-band absolute 
magnitude -29.5 < Mb < -27.0 and redshift 1.8 < z < 2.8 observed by MAMBO at 
250GHz, and a sample from Priddey et al. (2003a) of radio-quiet quasars with —29.2 < 



^Lr = lyL '■. L ° is the luminosity density at rest frame 4400A. we calculate L ° and Lr 

from the AB magnitude at rest frame 1450A in the discovering papers (Fan et al. 2000, 2001, 2003, 2004, 
2006a), assuming an optical spectral index of -0.5. 
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Mb < —27.5 and 1.5 < z < 3.0. Our higher redshift group (3.6 < z < 5.0) is drawn from 
the the Palomar Sky Survey (PSS) selected sample of Omont et al. (2001) (Mb < —27.0; 
3.9 < z < 4.6) and the SDSS sample of Carilli et al. (2001) (-28.8 < Mb < -26.1; 
3.6 < z < 5.0), both observed with MAMBO at 250 GHz. The 3o- detection limits of the 
MAMBO observations by Omont et al. (2001) and (2003) are ~1.5 to 4 mJy, and ~1.4 mJy 
by Carilli et al. (2001). The typical 3a detection limit of SCUBA observations reported by 
Priddey et al. (2003a) is ~ 10 mJy at 350 GHz, which corresponds to an upper limit of ~4 
mJy at 250 GHz (assuming a warm dust SED with temperature T^ = 47K and emissivity 
index P — 1.6). 

A total of 13 z~6 SDSS quasars were observed in the course of the work presented here, 
all with the VLA at 1.4GHz and eleven (all but J000552.34-000655.8 and J104845.05+463718.3) 
with MAMBO-2 at 250GHz (See Table 1 for the source hst). J104845.05+463718.3 was a 
pubhshed detection by MAMBO (Bertoldi et al. 2003a) and was observed by the VLA with 
a signal at 2.2a (Carilli et al. 2004). We re-observed it with the VLA, and combined all the 
data. 

MAMBO-2 at the IRAM 30m telescope is a 117-channel bolometer array. The Half 
Power Beamwidth (HPBW) of each pixel is 11" and the spacing between horns is about 20". 

The effective sensitivity is about 40 mJys^^^. The new observations were made in the winters 
of 2004-2005 and 2005-2006 during pooled observations using the on-off mode, wobbling by 
32" — 46" in azimuth and at a rate of 2 Hz. The target sources arc positioned on the most 
sensitive bolometer and the correlated sky noise is determined from the other bolometers 
and subtracted from the source bolometer. The median rms is ~0.8 mJy (Table 1). We 
reduced the data with the MOPSIC software package (Zylka 1998) and standard scripts for 
on-off observation data. 

The VLA observations were made with the A array, of which the maximum baseline is 
30 km. The sources were observed at 1.4 GHz with two Intermediate Frequency bands "IFs" 
and 50MHz bandwidth per IF. The corresponding Full Width at Half Maximum (FWHM) 
resolution is about 1.4". We observed each source for two to four hours to a typical rms noise 
level of ~16 //Jy beam~^ (Table 1). The data were reduced and images were made using the 
standard VLA wide field data reduction software AIPS. 



3. Results 

We present our MAMBO-2 results of 11 sources and the VLA results of 13 sources in 
Table 1. The optical properties are taken from Fan et al. (2006b) and presented in Col. (1): 
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the SDSS name, Col. (2): the redshift, and Col. (3): the AB magnitude at 1450A. Col. (4) 
lists the radio 1.4GHz surface brightness at the optical position. Col. (5), (6), and (7) list 
the nearest radio peak position and peak surface brightness for detected sources. Col. (8) 
presents the 250GHz flux densities in mJy. We mark detections in bold face. For the source 
J000552. 34-000655. 8, there is a strong (~Jy) radio source in the field which precludes deep 
VLA radio imaging. The rms on the 1.4GHz map is 130 /iJy, which is an order of magnitude 
higher than the others, yielding an extremely high upper limit of ~390 /xJy. Thus, we exclude 
this source in all of our analyses in this paper. 

We show the radio 1.4GHz continuum images of the 12 remaining sources in Figure 
1. We searched for > 3a peaks within a 0".6 radius from the optical quasar position in 
the radio map. This is a combination of the positional uncertainty of 0".3 between the 
radio and optical reference frames (Deutsch 1999) and the radio observation uncertainty of 
o'pos ~ ~ 0".eE. Six sources are detected at the > 3a level. According to the 1.4 

GHz radio source counts (Fomalont et al. 2006, in prep.), we expect 0.003 detections with 
Si.4GHz > 50 /iJy by chance, within the total search area around the 12 quasars. 

We obtain four detections among the 11 sources with MAMBO. Three of the detections, 
J084035.09+562419.9 J092721. 82+200123.7 and J133550.81+353315.8, are detected at > Aa. 
The fourth source, J081827.40+172251.8, is marginally detected at the ~ 3a level. As a 
comparison, the cumulative number counts of 250 GHz source is about 400 deg~^ with flux 
density S250GHZ > 2.4 mJy (Greve et al. 2004; Voss et al. 2006) which is the typical 3a limit 
of our observations. Thus we should expect 0.03 detections in our 11 fields by chance, given 
the 11" beam size of the MAMBO bolometer elements. 

We summarize the (sub)mm and radio observations of the six previously observed 
quasars in the z~6 SDSS sample (Petric et al. 2003; Priddey et al. 2003b; Bertoldi et 
al. 2003a; Robson et al. 2004; Carilli et al. 2004) in Table 2. Col. (1), (2), and (3) give the 
name, redshift, and 1450A AB magnitude, and Col. (4), (5), (6), and (7) the flux densities 
at 250 GHz, 350 GHz, 667 GHz and 1.4 GHz. 

Combining Tables 1 and 2, there are seventeen z~6 quasars that have been observed 
with MAMBO. Together with the results from SCUBA at 850 fim, there are eighteen z~6 
SDSS quasars that have (sub)mm observations at the 1 mJy sensitivity level and 8 of them 
are detected at > 3a. The detection rate is 44 ± 16%. The detection rate of optically bright 
quasar samples at redshifts 2 and 4 is ~ 30% in Omont et al. (2001, 2003) with an rms 
level of > 1.5mJy, and ~ 39% in Carilli et al. (2001) with a lower rms level of ~ lAmJy. 
Thus our detection rate at z~6 is slightly higher compared to these observations at lower 



^We require a signal to noise ratio (SNR) of 3 for detections here. 
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redshifts, but is consistent within the errors. 

Seventeen sources in the z~6 SDSS quasar sample have deep observations with the VLA 
at 1.4 GHz and 8 are detected at the > 3a level. J083643.85+005453.3 is the only >1 mJy 
radio source at 1.4 GHz among them. It was detected by Petric et al. (2003), as well as in the 
FIRST survey (White et al. 1997), leading to a weighted average flux density of 1.74±0.04 
mJ}!^!. Another source, J163033. 90+401209. 6, has not been observed with the VLA yet, and 
the FIRST catalog yields an upper limit of 0.44 mJy at the optical position (Bertoldi et al. 
2003). We thus exclude this source in the analysis of radio properties in the next section. 



To estimate the FIR luminosity of our z~6 quasars, we model the FIR continuum with 
an optically thin grey body (Eq. (1)), adopting a dust temperature of Td = 47K and emis- 
sivity index, j3 = 1.6, which is derived from the mean SED of high-redshift quasars (Beelen 
et al. 2006). We normalize the SED to the (sub)mm data and integrate the model SED 
from rest frame 42.5 to 122.5 fim to get the total FIR fiux0: 



This estimation of FIR emission is sensitive to the assumed dust temperature (Td). A higher 
temperature, eg. Td = 55K, will increase the estimated FIR flux by a factor of ~1.5. For 
non-detections (< 3a), we adopt as upper limits the larger value of either (a) the 2a rms or 
(b) the measured value at the optical position plus la rms for all the following calculations 
and plots. The FIR luminosities (Lfir) are given in table 4; most of the (sub)mm detected 
sources are very luminous in the FIR band, with implied Lfir ~ 10^^ (see Table 4). We 
estimated the FIR luminosities for all the low redshift comparison samples in the same way. 

We plot the FIR luminosity (Lfir) versus redshift in Figure 2, including the z~6 SDSS 
sample and all the comparison samples. In Figure 3, we present the histograms of the 
Lfir and Lfir/Lb distributions separately for different redshift intervals. For the (sub)mm 



3The 1.4 GHz flux density for J083643. 85+005453.3 is 1750 ± 40 /iJy in Petric et al. (2003) and 1530 ± 
150 /iJy in the FIRST catalog. 

"For J114816. 64+525150.3, we directly adopt the fitting results from Beelen et al. (2006) with Td = 55K 
and (3 — 1.6. 



4. Analysis 



4.1. 



The FIR emission 




(1) 
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detections at all redshifts, the Lpm values lie between 10^^-^ and 10^^-^ L0. The 1.5 < 
z < 3 group covers an even narrower range with most of the detected sources distributed at 
LpiR ^ 10^^ Lq These ranges are partially a result of the detection thresholds, i.e. a detection 
limit of l.SmJy at 250GHz corresponds to FIR luminosities ~ 10^^'^ for redshifts from 4 
to 6 and ~ 10^^ Lq at z~2. However, the Lfir/Lb distributions are similar for the (sub)mm 
detections at all redshifts (See Fig. 3 (b)). A similar conclusion was obtained by Omont et 
al. (2003) and Beelen et al. (2004), but we are extending it to the quasar sample at z~6. 
This fact suggests the optical to FIR spectra energy distributions (SEDs) of the luminous 
quasars in both optical and FIR bands do not evolve much from z~2 to 6. 

LpiR versus Lb is plotted in Figure 4, for quasars at all redshifts. Wc present correlation 
tests for (i) the z~6 sample only (z~6), (ii) samples at all redshifts (Total), and (iii) (sub)mm 
detections in all samples (Detection). We apply the General Kendall's tau test (Isobe et al. 
1986), taking Lb as the independent variable and Lpm as the dependent variable. The results 
are listed in Table 3 with Col. (1) the sample used. Col. (2) the sample size, and Col. (3) 
the Kendall's tau value. The probabilities (Pnuu) that no correlation exists between Lb and 
LpiR (the null hypothesis) are given in Col. (4). We take PnuU = 5% as the significance level 
below which the null hypothesis can be rejected. 

There is no correlation between Lb and LpiR for the z~6 sample only (PnuU ~ 44%). 
However, the PnuU value is 6% when all the samples are included, which suggests a marginal 
correlation. These results indicates that the FIR and optical emission of these (sub) mm 
observed quasars are correlated, but the scatter is large enough that the correlation is not 
seen over a narrow luminosity range. This correlation is even more significant when we do 
the test with only the (sub)mm detections at all redshifts (PnuU ~ 0.02%). But one should be 
careful with this result as there may be a number of observational or selection bialses in the 
sample of detections. For example, given the similar observational sensitivity level, the FIR 
luminosity threshold for the (sub)mm detections is decreasing with redshifts at z > 2 (see 
Figure 1). We apply linear regression to all of the (sub)mm observed samples (Total) and 
the subsample of detections (Detection), using the Expectation-Maximization (EM) method 
which can deal with data including upper limits (Isobe et al. 1986). The best fitting results 
are: 



Total : log (^^) = (0.21 ± 0.15)log(^) + (9.67 ± 2.03) 




(2) 



Detection : log (^^) = (0.40 ± 0.09)log(^^) + (7.82 ± 1.20) 




(3) 



The regression results show a non-linear relationship between the FIR emission and the 
quasars' optical emission with slopes much smaller than unity. We plot the fitting results 



-lo- 



in Figure 4. As a comparison, We calculate the typical FIR-to-B band luminosity ratios 
for local optical quasars based on the radio quiet and radio loud quasar templates in Elvis 
et al. (1994). The FIR luminosities are integrated directly from the template SEDs. The 
derived FIR-to-B band luminosity ratios (Lfir/Lb) are 0.29 and 0.38 for the radio quiet 
and radio loud templates respectively, and are ploted as Lfir = 0.29Lb (solid line) and 
LpiR = 0.38Lb (long dashed line) in Figure 4. According to Figure 4, most of the (sub)mm 
detections have larger Lfir/Lb value than that of the local optical quasar templates. For 
the eight (sub)mm detected quasars at z~6, six of them have FIR emission stronger than 
that predicted from the radio quiet template (by factors from ~ 1.5 to > 5), and only one 
source J081827.40+172251.8 is consistent with the template within errors. 



4.2. The Radio emission 

The rest frame 1.4GHz luminosities (LRaJ^ for the 17 quasars observed with the VLA 
in the z~6 sample are calculated, assuming a power law (f^, ~ i/") radio SED and a radio 
spectral index of a = —0.75 (Condon 1992). Since all the sources are unresolved on the radio 
maps (Figure 2), we adopt the peak surface brightness in Table 1 as the 1.4GHz flux density 
for our VLA detections and take the larger value of either the 2a rms on the radio map or 
the measured value at the optical position plus Icr rms as the upper limits for non-detections. 
The rest frame radio to B band optical luminosity density ratios R^ 4 = L,^ i,4ghz/L ^^^q^ 
(Cirasuolo et al. 2003; 2006) are calculated and listed in Table 4. L -is the luminosity 

^ ' ' v,4400A 

density at rest frame 4400A (See the footnote in Section 2). We adopt R = Li,^5ghz/L ^^q^^ = 
10 (Kellermann et al. 1989) to separate radio loud and radio quiet sources. This corresponds 
to R* 4 ~ 30 by assuming the above radio spectral index and converting the rest frame 5 
GHz flux density to 1.4 GHz. 

We plot the rest frame 1.4 GHz radio luminosity (LRad) vs. B band luminosity (Lr) 
in Figure 5. The typical 3a detection limit of the VLA observation is ~ 50 /iJy (dotted 
line), which is roughly 10 times deeper than the FIRST survey (solid line). Thus our VLA 
observations are sensitive enough to detect any radio loud source with B band luminosity 
> 10^^ L0 at z~6 for all of the 17 sources. We plot the R^ 4 distribution of the 17 z~6 
quasars in Figure 6. None of the 13 sources with new radio data reported in this paper 
is radio-loud. For the whole sample of 17 z~6 quasars, there is one marginally radio-loud 
object, J083643.85+005453.3, with R* 4 ~ 40. 



■''LRad = i^Li/,i.4GHz in rest frame 



The deep VLA observations also enable us to investigate the correlation between radio 
and FIR emission of the z~6 quasars. Yun et al. (2001) employ a q parameter to quantify 
the ratio of the FIR and radio luminosities: 



q = log( 



FIR 




fl.4GHz 



)■ 



(4) 



3.75 X 10^2 Wm-: 



Wm-2 Hz-i 



Yun et al. (2001) studied the infrared-selected galaxies in the IRAS 2Jy sample and found 
a q value of 2.34 for typical star forming galaxies. We calculate the q values for the 12 z~6 
quasars that are detected at either radio or (sub)mm wavelengths. 

We also plot Lpm vs. the rest frame 1.4GHz radio luminosity density (L,^^i.4ghz) in 
Figure 6, comparing the z~6 quasar sample to typical star forming galaxies in the local 
universe (Yun et al. 2001). The solid line in Figure 6 corresponds to the typical q value of 
2.34 in star forming galaxies, and the dotted lines denote excesses that are 5 times above 
and below this typical correlation (Yun et al. 2001). 

There are four sources in the z~6 quasar sample that are detected in both the radio 
and (sub)mm regimes. One of them, J081827.40-I-172251.8, has a small q value that falling 

beyond the range defined by typical star forming galaxies, indicating the dominance of AGN 
power in the FIR to radio SED of this source. However, the FIR-to-radio ratios of the other 
three sources, though slightly above the mean value of q = 2.34, are consistent with star 
forming galaxies (Figure 6 and Table 4). 



The heating sources of the FIR emitting dust in quasars are studied and the contri- 
butions from host galaxy star formation arc discussed in a number of papers (eg. Haas et 
al. 2003; Omont et al. 2003; Schweitzer ct al. 2006). Star formation is believed to be 
the dominated dust heating resource in the local infrared luminous quasars that are located 
in ULIRGs (Hao et al. 2005). Schweitzer et al. (2006) studied the connection between 
FIR emission and PAH/low excitation fine-structure emission lines of the local PG quasars. 
Their results suggest that star formation is responsible for at least 30% to the average lo- 
cal quasar FIR emission and that this contribution tends to increase with FIR luminosity. 
The z~6 SDSS quasars we discussed in this paper belong to the population of the brightest 
quasars in the optical. About 1/3 of them are detected at (sub)mm wavelength. According 
to Figure 4, most of these detections have FIR luminosities stronger than the predictions 
from local quasar templates, indicating extra emission from warm dust in these sources com- 
pared to the typical optical quasar emission. Moreover, the relationship between the FIR 



5. 



Discussion 
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and B band optical emission is non-linear and significantly scattered, only manifesting it- 
self when a larger luminosity range is considered (See Section 4.1). Bcelen (2004) studied 
a number of optically bright quasar samples at z > 1.5 that had (sub)mm observations, 
and first report this weak correlation between FIR and B band optical emission, namely 
logLpiR = (8.36 ± 0.90) + (0.33 ± 0.07) x logLe (See also Cox et al. 2005). The resuhs are 
consistent when we include our new observations at z~6. Additionally, most of the (sub)mm 
detections at z~6 have FIR-to-radio ratios or upper limits within the range occupied by star 
forming galaxies (See Table 4 and Figure 6). 

One possible explanation for all of these facts is that star formation is happening on 
a massive scale in the host galaxies of these FIR luminous z~6 quasars and dominating 
the heating process of the FIR emitting warm dust. The implied star formation rate is 
~ lO^M0yr~^. This idea is supported by the detections of CO in the high-redshift quasars 
that have strong FIR emission (Lpm > 10^^ L©). CO hue emission was detected in several 
FIR luminous quasars at z>4, including the highest redshift quasar J114816. 64+525150. 2 
(eg. Riechers et al. 2006; Bertoldi et al. 2003b; Walter et al 2003,2004). The CO detections 
indicate huge amount of molecular gas in the host galaxies of these FIR luminous quasars 
which are the requisite fuel for star formation. Moreover, Riechers ct al. (2006) found that 
the FIR-to-CO relation of the high-redshift FIR luminous quasars are consistent with that 
defined by ULIRGs, starforming galaxies and submillimeter galaxies. 

One may argue that the strong FIR emission can be also processed by the AGN power 
since the dust torus geometry of these sources arc still unkown. This possibility cannot be rule 
out given the limited data we have for these sources, but detailed dust models are required 
to explain the results we hst above. One fact is that the hot dust emission of the SDSS 
z~6 quasars probed by recent Spitzer observation are similar to those of the local quasar 
templates as we mentioned in Section 1 (Jiang et al. 2006a). Four (sub)mm detections are 
included in this Spitzer observed sample, and no different properties are found in the optical 
to near infrared SEDs of these sources. This may not support the idea of a quite different 
geometry of the dust torus. 

However, we should mention that there are still large uncertainties in the estimations 
of the FIR luminosities for these (sub)mm detected z~6 quasars. The FIR luminosities for 
some of the sources can be overestimated given the poor data at (sub)mm wavelengths. Thus 
further studies on these (sub)mm detected z~6 quasars are required, including observations 
of the submm continuum at higher frequencies and emission lines of CO and other molecules, 
such as HCN, as well as the C and C^ fine structure lines. 

Another interesting question is how the current radio observations constrain the radio 
loud fraction (RLF) at z~6. The RLF of optically selected quasar samples is ~ 10%, as 
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quoted in many papers (Kellermann et al. 1989; Ivezic et al. 2002, 2004; Cirasuolo et al. 
2003, 2006). Jiang et al. (2006b) studied the FIRST data (Becker et al. 1995) for SDSS 
quasars with redshifts up to 5, and suggested that the quasar RLF increases with optical 
luminosity and decreases with redshift, namely log ((RLF) /(I — RLF)) = (— 0.112±0.109) + 
(-2. 196 ± 0.269)log(l + z) + (-0.203 ± 0.026) (M2500 + 26). This result gives an RLF of ~ 37% 
at z = 0.5, ~ 11% at z = 2, and ~ 2% at z = 6, given the rest frame 2500A absolute AB 
magnitude M2500 = —27.3 — the typical value of the z~6 SDSS quasar sampl^. 

There is one marginally radio loud source in this z~6 sample of 17 sources. This result 
argues against a high RLF, i.e. > 20%, for the current SDSS quasar sample at z~6. Thus 
it is roughly consistent with the result of Jiang et al. (2006b). However, the sample size is 
still too small to set a strict constrain, i.e. differentiate an RLF between ~ 2% and ~ 10%. 
Thus to provide a better test of the redshift evolution of the RLF; we would need a sample 
three times larger to usefully constrain the RLF at z~6. 

We should also mention the recent discovery of a radio loud quasar at z>6, FIRST 
J1427385+331241(McGreer et al. 2006), with a FIRST flux density of 1.73mJy. This source 
is the most radio luminous quasar known at z~6, with RJ^ > 100. However, it cannot be 
included in our radio loud analysis at z~6, as the selection criteria of this source are quite 
different from that of the SDSS quasars and its UV/Optical emission is just beyond the 
SDSS detection limits (McGreer et al. 2006). 



6. Conclusion 

In this paper, we present new results of millimeter and radio observations of a sample of 
z~6 quasars. These quasars are selected from the SDSS survey and observed with MAMBO- 
2 at 250 GHz and the VLA at 1.4 GHz. We obtained three > 4cr detections and one ~ 3cr 
detection out of 11 observed sources by MAMBO-2 and six radio detections out of 13 observed 
sources by the VLA. 

We combine our new millimeter and radio results of the SDSS z~6 quasars with results 
from the literature and discuss the FIR and radio properties of the optically selected quasars 
at z~6. Our conclusions are as follows. 

• Eight out of 18 z~6 optically selected quasars are detected in the (sub)mm regime. 
This indicates a (sub)mm detection rate of 44% at z~6 at mJy sensitivity. Within 



®Thc M2500 is estimated from the absolute AB magnitude at rest frame 1450A (M1450, Fan et al. 2006b). 
For the current sample of 19 z~6 SDSS quasars with M1450 from -26.2 to -27.9, we adopt a typical value of 
M1450 — —27.0 and an optical spectral power law index of -0.5. 
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the errors, this is consistent with the ~30% (sub)mm detection rate at lower reshift 
(eg. CariUi et al. 2001; Omont et al. 2001; 2003). The observational data imply FIR 
luminosities ~ 10^^ Lq in the (sub)mm detected sources. 

• We compare the distribution of FIR luminosities and FIR to optical ratios between the 

z~6 SDSS sample and (sub)mm observed optically bright quasars (Lb > 10^^'^ Lq) at 
lower reshifts. The distributions of the FIR-optical ratio are similar for different red- 
shift groups, which suggests that the average optical-to-FIR SED of optically bright 

quasars is independent of redshift. 

• We extend the quasar FIR-to-optical correlation study to the z~6 SDSS sample. No 

correlation is found with the z~6 sample only. However, a correlation (albeit with large 
scatter) can be seen when all the samples extending from z=1.5 to 6.42 are included. 

• We also discussed the FIR-to-radio ratios of the z~6 quasars, by comparing them to 

the typical correlation defined by star forming galaxies. Three of the four sources that 
are detected in both the millimeter and radio bands have FIR-to-radio ratios within 

the range defined by star forming galaxies. 

• We found no strong radio sources with > 30 among the new SDSS sources of 

bright z~6 quasars. In the whole z~6 sample, only one radio detection has a radio to 
optical ratio Rl^ ~ 40 and no source has R^^ > 100. These data are consistent with, 
although do not set strong constraints on, the recent conclusion of a decreasing radio 
loud quasar fraction with increasing redshift (Jiang et al. 2006b). 

These results give a view of the general FIR through radio properties of the z~6 SDSS 
quasars. The data are consistent with the idea that massive starbursts may exist in the 
host galaxies of the strong (sub) mm detections at z~6 and contribute to the FIR and radio 
emission. These strong (sub)mm sources provide the only candidates to search for CO and 
C"^ into the epoch of reionization, and to test the idea of co-cval SMBH and host galaxy 
formation. We may expect to go an order of magnitude deeper in a few years time with 
the coming instruments of the Expanded Very Large Array (EVLA) and the Atacama Large 
Milhmeter Array (ALMA). 
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Fig. 1. — VLA images of 12 z~6 quasars at 1.4GHz at 1.4" resolution (FWHM). The param- 
eter 'con' denotes the value for the contour increment in each map in units of /^Jy/Beam. 
The contour levels are (-3, -2, 2, 3, 4, 5, 6, 7, 8, 9, and 10) xcon. The eUipse indicates the 
beam in each case, and the cross marks the optical position. 
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Fig. 1 - Continued 
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Fig. 2. — The logarithm of the FIR luminosity (logLFm) versus redshift for different samples. 
The filled symbols represent detections with la errors. The open symbols with arrows denote 
upper limits: we adopt the larger value of either the 2a rms or the measured value at the 
optical position plus la rms as upper limits for non-detected sources. The dashed and dotted 
lines represent the typical 3a detection limits of MAMBO at 250 GHz and SCUBA at 350 
GHz, respectively, namely S250GHZ = 2.4 mJy (this work) and Sssoghz = lOmJy (Priddey et 
al. 2003a). 
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Fig. 3. — (a) The distribution of the logarithm of the FIR luminosity for quasars at different 
redshifts. The z~6 quasars are plotted as black lines. The 1.5 < z < 3 group (red lines) is 
combined with the samples from Omont et al. (2003) and Priddey et al. (2003a), while the 
3.6 < z < 5 group (blue lines) is from Omont et al. (2001) and Carilli et al. (2001). For 
all the samples, the shaded areas represent detections and dotted lines denote upper limits: 
the upper limits are calculated as in Figure 2. (b) The distribution of the logarithm of the 
FIR-optical ratio (logLFiR/Le) for quasars at different redshifts. 
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Fig. 4. — The logarithm of the FIR luminosity is plotted against the logarithm of the B band 
luminosity for all the (sub)mm observed quasars. The symbols are the same as in Figure 
2. The dashed line is the linear regression result for the (sub)mm observed samples at all 
redshifts including upper limits, while the dotted line is the result using the detected objects 
only. The solid and long dashed lines represent the FIR-to-B band luminosity ratios of local 
radio quiet and radio loud quasar templates in Elvis et al. (1994), respectively. 
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Fig. 5. — The logarithm of the rest frame 1.4GHz radio luminosity (logLRad) is plotted 
against the logarithm of the B band luminosity (logLe) for the z~6 quasars observed with 
the VLA at 1.4GHz. The filled squares represent detections with Icr errors, and the open 
squares with arrows denote upper limits: the upper limits are calculated as in Figure 2. The 
dashed hues represent rest frame radio-to-optical ratios (R1.4) of 100, 30 (separation of radio 
loud and quiet), and 1. The dotted line denotes the typical 3cr detection limits of our VLA 
observations, while the solid line shows the Scr detection hmit of the FIRST survey. 



-25- 




11.0 11.5 12.0 12.5 13.0 13.5 

log Lpir(Lq) 



Fig. 6. — The logarithm of the rest frame 1.4 GHz radio luminosity density (logL,^^i.4GHz) 
versus the logarithm of the FIR luminosity (logLpm) . The filled squares represent the de- 
tections in both radio and (sub) mm bands, and the open squares with arrows represent the 
upper limits of the non-detections in either radio or (sub) mm. The upper limits are cal- 
culated as in Figure 2. The crosses represent the IRAS 2Jy sample of galaxies in Yun et 
al. (2001) and the dashed line indicates the typical radio-to-FIR correlation in star forming 
galaxies, with correlation parameter q = 2.34 (Yun et al. 2001). The dotted lines represent 
excesses five times above and below the typical q value. 
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Table 1: 1.4GHz and 250GHz results of 13 z~6 quasars. 



(1) 


(2) 


°^14B0A 

(3) 




radio detection at 1.4 GHz 




Q 

'~J250GHz 

mjy 

(8) 


o b 
Sl.4opt 

TO —1 

(4) 


h m s 

(5) 


o ' " 

(6) 


Q C 

'-'peak 
TO —1 

(7) 


J000239. 39+255034.8 


5.80 


19.02 


81+14 


00:02:39.36 


25:50:34.9 


89'^ 


0.20+0.88 


J000552. 34-000655.8 


5.85 


20.03 


40+130 










J081827.40+172251.8 


6.00 


19.34 


119±12 


08:18:27.39 


17:22:51.8 


123 


1.45+0.49 


J084035.09+562419.9 


5.85 


20.04 


12±9 








3.20+0.64 


J092721.82+200123.7 


5.79 


19.87 


33±14 


09:27:21.82 


20:01:23.1 


45 


4.98+0.75 


J104845.05+463718.3 


6.20 


19.25 


6±11 








3.00+0.40* 


J113717.73+3549.56.9 


6.01 


19.63 


9±17 








0.10+1.13 


J125051. 93+31.3021. 9 


6.13 


19.64 


37+21 








0.07+0.90 


J1.33550. 81+35.3315.8 


5.95 


19.89 


35+10 


13:35:50.81 


35:.33:15.9 


35 


2.34+0.50 


J141111.29+121737.4 


5.93 


19.97 


44±16 


14:11:11.29 


12:17:38.0 


61 


1.00+0.62 


J143611.74+500706.9 


5.83 


20.16 


6±16 








-0.21+1.14 


J160253.98+422824.9 


6.07 


19.86 


53±15 


16:02:53.95 


42:28:24.9 


60 


1.82+0.86 


J162331.81+311200.5 


6.22 


20.13 


24±31 








0.17+0.80 



■^Result taken from Bcrtoldi ct al. (2003a). 

''Surface brightness at the optical quasar position plus the la rms on the radio map. 
'^Position and surface brightness of the detected radio peak. 
''The detections are marked as boldface. 



Table 2: P ublished results from the literature. 



SDSS name 

(1) 


z 

(2) 


m 

14,-)0A 

(3) 


S250GHZ 
mJy 
(4) 


S350GHZ 
mJy 
(5) 


S667GHZ 
mJy 
(6) 


S 1.4GHz 

ujy 

(7) 


J083643.85+005453.3 


5.82 


18.81 


-0.4+1.0=' 


1.7+1.5*^ 


-24+10'' 


1740+40f'9 


J103027.10+052455.0 


6.28 


19.66 


-1.1+1.1* 


1.3+1.0'' 


-21+10'' 


-3+20* 


J104433.04-012502.2 


5.74 


19.21 




6.1+1.2'' 




-15+24* 


J114816.64+525150.3 


6.42 


19.03 


5+0.6= 


7.8+0.7^^ 


24.7+7.4^ 


55+12"^ 


J130608.26+035626.3 


5.99 


19.55 


-1.0+1.0* 


3.7+1.0'' 


-7+14'' 


14+21* 


J163033.90+401209.6 


6.05 


20.64 


0.8+0.6'^ 


2.7+1.9'' 


15.4+9.6'' 


< 440= 



"Petric et al. (2003). 
''Priddey et al. (2003b). 
'^Bertoldi et al. (2003a). 
''Robson et al. (2004). 
'^CarilH et al. (2004). 

•^weighted average of the results in Potric et al. (2003) and the FIRST survey. 
^The detections are marked as boldface. 
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Table 3: Correlation tests 



Group 


Number 


Kendall's tau 


Pnull 


(1) 


(2) 


(3) 


(4) 


z~6 


18 


0.2222 


0.4402 


Total 


208 


0.1246 


0.0575 


Detection 


64 


0.6310 


0.0002 



Table 4: Derived parameters of the z^6 quasars 



SDSS name 


logLs 


logLpiR 


logLRad 




q 




(Lq) 


(L©) 


(L©) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


J000239.39+255034.8 


13.20 


<12.62 


7.89 


2.40 


<1.31 


J081827.40+172251.8 


13.09 


12.53 


8.06 


4.55 


1.04 


J083643. 85+005453.3 


13.28 


<12.66 


9.18 


38.68 


<0.05 


J084035. 09+562419.9 


12.79 


12.88 


<7.27 


<1.46 


>2.18 


J092721.82+200123.7 


12.86 


13.08 


7.59 


2.65 


2.06 


J103027. 10+052455.0 


12.99 


<12.71 


<7.62 


<2.05 




J104433. 04-012502. 2 


13.11 


12.83 


<7.62 


<1.55 


>1.79 


J104845. 05+463718.3 


13.15 


12.84 


<7.42 


<0.90 


>1.99 


J113717.73+354956.9 


12.98 


<12.72 


<7.50 


<1.65 




J114816.64+525150.2 


13.26 


13.18 


7.77 


1.60 


1.98 


J125051.93+313021.9 


12.98 


<12.62 


<7.75 


<2.87 




J130608.26+035626.3 


13.00 


12.61 


<7.60 


<1.91 


>1.59 


J133550.81+353315.8 


12.86 


12.74 


7.51 


2.14 


1.80 


J141111. 29+121737.4 


12.83 


<12.58 


7.75 


4.01 


<1.41 


J143611. 74+500706.9 


12.74 


<12.73 


<7.45 


<2.48 




J160253. 98+422824.9 


12.89 


<12.76 


7.76 


3.62 


<1.58 


J162331.81+311200.5 


12.80 


<12.56 


<7.80 


<4.86 




J163033.90+401209.6 


12.58 


<12.51 









Note. — Col. (1), the name of the source, Col. (2), the B band luminosity (see the footnote 
in Section 2 for the calculation), Col. (3), the rest frame FIR luminosity. Col. (4), the rest 
frame 1.4 GHz radio luminosity. Col. (5), the radio-optical ratio R14, and Col. (6) the 
FIR-to-radio correlation parameter q. 



